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Abstract — Microgrid is an aggregation of distributed 
generators (DGs) and energy storage systems (ESS) through 
corresponding power interface, such as synchronous generators, 
asynchronous generators and power electronic devices. Without 
the support from the public grid, the control and management of 
an autonomous microgrid is more complex due to its poor 
equivalent system inertia. To investigate microgrid dynamic 
stability, a small -signal model of a typical microgrid containing 
asynchronous gener-ator based wind turbine, synchronous 
diesel generator, power electronic based energy storage and 
power network is proposed in this paper. The small-signal model 
of each of the subsystem is established respectively and then the 
global model is set up in a global reference axil frame. 
Eigenvalues distributions of the microgrid system under certain 
steady operating status are identified to indicate the damping of 
the oscillatory terms and its effect on system stability margin. 
Eigenvalues loci analysis is also presented which helps 
identifying the relationship among the dynamic stability, system 
configuration and operation status, such as the variation of 
intermittent generations and ESS with different control 
strategies. The results obtained from the model and eigenvalues 
analysis are verified through simulations and experiments on a 
study microgrid system. 

Index Terms —Distributed generators, dynamic stability, 
energy storage system, microgrid, small-signal model. 

I. INTRODUCTION 

Microgrid can be considered as an advanced organi-zation 
of distributed generation systems with less infl u-ence from 
the stochastic distributed generations (DGs), such as PV 
arrays and wind turbines, through proper control of energy 
storage systems (ESS) and dispatched loads [1], [2]. 

Microgrid plays important roles in the effective utilization 
of renewable energies and the stable operation of the public 
power grid [3], [4]. It makes the massive development of the 
renewable ener-gies possible. 

Normally, DGs can be connected into microgrid through 
power electronic converters or rotating machines [5]. The 
former responds rapidly to the power control and make 
mi-crogrid operation more flexible; however, they are 
potentially susceptible to the oscillation under system 
disturbances [6]. 

While the latter usually involves synchronous generators 
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equipped with excitation and governor systems, and 

asyn-chronous wind turbine (AWT) of small or medium 
capacity. 

While the latter usually involves synchronous generators 
equipped with excitation and governor systems, and 

asyn-chronous wind turbine (AWT) of small or medium 
capacity. 

A typical microgrid usually has two kinds of operation 
modes, i.e., autonomous operation mode and grid -connected 
operation mode. To guarantee the reliable power supply for 
essential loads, the smooth transition between these two 
modes is usually required [7], [8]. Without the support from 
the main grid, the control and management of the autonomous 
operation microgrid is often more complex than the 
grid-connected one. Since the equivalent physical inertia of 
microgrid is quite smaller compared with the main grid [4], 
some medium or even small interruptions such as the output 
power fluctuation of PV arrays or wind turbines, and the cut in 
or cut off of local loads may result in power quality and 
stability issues. Specifically, the frequency and voltage 
quality would decline, and system stability would deteriorate. 

With the development of power electronics, energy storage 
can be controlled to release or absorb active/reactive power 
flexibly. When ESS is integrated into microgrid, it can 
main-tain system instantaneous power balance and improve 
the dy-namic performance through proper energy 
management strate-gies, which equals to system equivalent 
inertia enhancement [9]. Therefore, the impact of DGs’ power 
fluctuations and other in-terruptions on system stability and 
power quality can be miti-gated greatly, which is of vital 
importance for microgrid in au-tonomous operation mode. 
Many energy storage technologies and their applications in 
microgrids were studied recently in-cluding Li-ion batteries, 
flywheel energy storage (FES), super-capacitors, 
superconducting magnetic energy storage (SMES) and 
vanadium redox battery (VRB), etc. [10], [11]. 

One of the important concerns in the reliable operation of 
mi-crogrid is small-signal stability [8]. In conventional power 
sys-tems, stability analysis is well established with standard 
models of synchronous machines, governors and excitation 
systems of varying orders that are known to capture the 
important modes for particular classes of problems. This does 
not yet exist for microgrids and may be hard to achieve 
because of the wide range of power technologies that might be 
deployed. Therefore, full- order models of microgrids with 
typical DGs are required. The small-signal models containing 
synchronous generators, in-verter based generations and 
power network were established in [5] and [12], and the 
dynamic characteristics in autonomous operation mode were 
analyzed based on the individual models. The global 
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small-signal model of the inverter based microgrid was setup 
in [8] to design the controller of inverter based DGs. 
However, as an important DG, asynchronous wind turbine 
was not included in the model. 

In this paper, a small- signal model of a typical microgrid is 
presented. It will contain the major types of power interface 
de- 


Mlcra-Grid 



Fig. 1. Single-line diagram of a typical microgrid. 


vices of DGs exist nowadays such as synchronous generator, 
asynchronous generator and inverter. Each of the small-signal 
models of the individual DGs will be established respectively 
together with that of the power network, and then the global 
model will be set up in a global reference axil frame. Once the 
small -signal model has been constructed, eigenvalues are 
iden-tified to indicate the damping of the oscillatory terms 
under dif-ferent operation conditions. The analytical nature of 
this exam-ination then allows further investigation of the 
relation between stability and system configurations, such as 
DGs capacity and system operating state. A sensitivity 
analysis is then conducted by tracking eigenvalues loci under 
certain conditions, which provides the difference between the 
situations with/without ESS on system stability. The 
theoretical analysis is verified with both simulations and 
experiments of a typical microgrid. 


the transfor-mation method shown in Fig. 2 and formulated in 

(l)[13j. 

As shown in Fig. 2, the (/-axis and tf-axis is taken as the 
global reference frame which is defined on the microgrid 
network ro-tating at angular frequency of uj . di -axis and ( h 
-axis, dz -axis and Qz -axis, d -axis and Q -axis are the local 
reference frames of DG1, DG2 and ESS rotating at uj , uj and 
uj , respectively. 



Fig. 2. Global and local rotating reference frames of the 
microgrid. 


bnin—wsc) denotes the angle between the local reference 
frame of dn— ( ln(n= 12a) and the global reference frame, 
respectively: 



In (1), \fnfdq It is defi ned as the vector of/, which repre-sents 
the voltage or current component in the nth local reference 
frame, [fqfg. Jg is the vector off in global reference frame, and 
l/^Jr is the steady-state value off in the nth local refer¬ 


The rest of the paper is organized as follow. Section II 
de-scribes the mathematical model of a typical microgrid 
system. Section III analyzes the small -signal dynamic 
stability by eigen-values loci tracking. Results and 
discussions of the simulation and experiment are reported in 
Section IV. Conclusion is stated in Section V. 


II. MATHEMATIC MODEL OF MICROGRID 

Fig. 1 shows the single-line diagram of a typical microgrid. 
It comprises, instead of local loads, DG1 representing an 
AWT with stall regulation, DG2 representing a synchronous 
gener-ator equipped with excitation and governor control 
systems, and an ESS based on a voltage-sourced converter 
(VSC) with ac-tive/reactive power controller. It is connected 
to a distribution network at the point of common coupling 
(PCC), and operates autonomously when the main grid is 
fault. 

The state-space model of DG1, DG2, ESS and power 
net-work are represented on their individual local reference 
frame, respectively. A common reference frame is chosen as 
the global reference frame of the microgrid system, and all of 
the subsys-tems are shift to the global reference frame using 


ence frame. 


A. Mathematic Model of DG1 


DG1 is mainly composed of wind turbine, transmission 
system, asynchronous generator and pitch control system. The 
electrical system of the three-phase symmetrical 
asynchronous generator in its arbitrary reference frame of 
di—fh is given by 

(/ 

\/ u ' = El w + E 1 — (l w ) (2) 

where 



are the voltage and current vectors of the stator windings 
(dllsdl an)d the rotor windings A/lrr/l, respectively, and 
matrices E and F can be derived in [13]. 

The rotor mechanical model of DG1 is represented as 


d/jj r 


211 


(F, - T) 



where ujb and uj are the base electrical angular velocity and 
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the rotor electrical angular velocity, respectively. H is the 
combined inertia constant of rotor and loads. T is the 
mechanical torque, and T is the electromagnetic torque given 
by e 

_ (4) 

where [ is the mutual inductance between stator and rotor 
windings. 


The characteristics of wind turbine are stated as 

{ ~ (5) 

where A is the ratio of rotor blade tip speed to wind speed. Re 
is the turbine swept radius. ujwmd is the angular velocity of the 
wind turbine. \ wind is wind speed. Tv j s the mechanical 
output torque of wind turbine. P is air density, f'r is the 
torque coeffi-cient. ft is the blade pitch angle [14]. 


Wind turbine is connected to the asynchronous generator 
through a coaxial shaft where the mechanical model can be 
regarded approximately as a first-order time-delay system due 
to its large physical inertia, as shown in (6): 


dT 

tit 




where Th i s the inertia time constant. 



The concept of slip s is defined in synchronous rotating 
ref-erence frame by setting uj=uj : 


where uk. is the synchronous electrical angular velocity. Then 
the model in synchronous rotating reference frame of dvJh is 
transformed into the global reference frame by (1). 


The mathematical model of DG1 can be set up by (2)-(7). 
The small- signal state space model of DG1 in the global 
reference frame is 

AT«’» = 4 u , A + *■' A»>" + *" A ( 8 ) 


A*"* = [A** tI A i* 4 A*,,, At,vi Aw, AT) t 
— i A\ /f AV* l r 

Art"' = |AV rf i av;„ Aiy r 
_ KtCrKKind ji; 

» It* — rp ^ ' ivtnd' 

• H 

B. Mathematic Model of DG2 

DG2 consists of a prime mover su^h as diesel engine and 
gas turbine, and a three-phase synchronous generator with 
excita-tion and governor control systems. The electrical 
system of the synchronous generator in rotor reference frame 
d. 2 — ( l 2 is 

(9) 


IT = or +H — (T) 

* L 


where 


V‘ = [V q2 V d2 V klq2 V k2q2 V kd2 V^ d2 ] T 
I* — [* 4 2 ii2 4-1*2 U-2,2 ikd‘2 ifd 2 ] 


are the voltage and current vectors of the stator windings 
(f 2 , ^ 2 ) ? the damper windings (k 1 q2 , k2q 2. k 4:2 ) and the field 
winding (fd2), respectively. The matrices G and H can be 
derived from the methods shown in [13]. 


The mechanical model of the rotor of DG2 is given by 



Fig. 3. Schematic diagram of ESS 


where uj is the rotor electrical angular velocity. H s is the 
com-bined inertia constant of rotor and loads. T M i s the 
mechanical torque. T rs i s the electromagnetic torque given by 


+ if 42 + 

- A" m 4 iif2 + 4 L #2 + ) ij2 (11) 


where Xmd and Xm. ^re the magnetizing inductances. 

The mathematical model of DG2 can be established by 
(9)—(11), and transformed to the global reference frame 
pro-vided that the synchronous generator operates at the 
grid-con-nected mode. It is convenient to linearize the state 
equations in the synchronous rotating reference frame with 
the relationship 
of _ and | given by 


dS m 
■ v 

~w~ 


^,( 12 ) Where u)b wh 


Combined(12) with(l) a the small-signal model ofDG2in 
the global reference frame are 


Xx** = + JJAtf 


03) 


Where 


A?.-* 


= A i* 2 AtAikiji Aik^2 Aikj'j Aijj% AJ n 1 
= [AV£AV* 2 ] T 
= |A V Id2 AT M \ T . 


C. Mathematic Model of ESS 


Fig. 3 illustrates the schematic diagram of ESS, which consist 
of battery bank and power electronic circuit with corresponding 
control system, and filters. The filters can be equivalent as the 
series connection of a resistor (r) and an inductance C L) for each 
phase where R denotes the lump resistor of power line and 
converter, and the effect of the dead-time of IGBT bridges. Xwc 
is the DC -bus voltage and can be considered constant usually 
for dynamic analysis, which explains the concept of VSC. 

In the local reference frame of d— ( J , the power circuit 
equations is given by 
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dt J 


a 

T ^ 

i t R 

UJi - T 


t'dc. 

J L'Vj 


1 

+ L 


V dc - 


L V, C 


t rN 

[ dc 

yN 

V f c 


where V’a and V’a are the voltage components of VSC in 
dc—flc axis frame. Vdc and Vc are the voltage components 



Fig. 4. External active/reactive power controller. 


of AC-bus in dcMlc axis frame, id: and ic are i\]p current 
components of VSC in dcMlc axis, respectively. 

Based on (14), the small-signal model of the power circuit 
of ESS is represented by 


Ax *9 = A c Ax e * -f 9 *Av c + + 9 * U AS C 


(15) 


Where 


AM* 

= |a,; c a«»j t 


Av c 

=[av; c av* c } 

T 

Au c 

= [av; av;] t 


K 

= v, e - v* - 


Vi 

= V* - V* + 

/ ^ 1 j • 


In the global reference frame,A*., can be derived from [12]. 

A novel control strategy is put forward to control the VSC 
of ESS. It is composed of three controllers including external 
active reactive power controller, inner current decoupling con¬ 
troller, and voltage feed-forward controller. The power con¬ 
troller generates the reference currents for the inner current con¬ 
troller according to the power requirements locally or remotely. 
Ihe inn^r current controller md voltage feed-forward controller 
th en gen erate th e r e ference v olta ge s (and \^) for VSC . 
As shown in Fig. 4, the phase lockloop (PLL) is used to 
recognize the instantaneous phase of the electric network. The 
mathematical model of the PLL is represented by 

Auj, = itAK, + —) = h',A.j K,Ai r (16) 

where K pa A" t are the PI coefficients of the PLL controller. 

The small-signal model of the control system of ESS in the 
global reference frame is given by 





+ 


fdn 

fqn 


T = 


cos S n 
sin 6 n 


— sin b n 
cos 6 n 



and_are the states of PI regulator of the current controller 

in_axis._and are the active and reactive power 

an 


output of ESS. and are the filtered values of d 


i ii i 


an 

II .||| and | are the filtered values of 

and | 

. -J d 

It — l D~ 

respectively. 


,(15) 

Based on the relationship between 

and 

can 

be rewritten as 

Ax'* = + JSAif + *;; k Aj'™ 

4- fA 

S' 

^=1 

r 

g 


D. Mathematic Model of Network 

As shown in Fig. 1, the model of network ina&cc stationary 
reference frame is represented by 


V 1 + 


+ 


Ihi li + 


L\ ii n 

it 

Lzii 13 


it 

dt 


__ j® 


- „ FCC ' 


l: FCC _ 


i" t* 3 



where/?,, and L a are the line impedance of the h th branch^ 1 
is the ^th branch current, v?™ ^ and are the voltage 

vectors ofPCC, Bus-I a Bus-2wdBus-3 a respectively £2, is the 
equivalent capacitor of the loads, and 



wh ere i* i s the current v ector o f DG2. 

The small-signal model of the power network of the micro¬ 
grid system in the global reference frame is represented by 


Ax rt ct = A „ rt A r t j r f H- 9“^ t A v + 9^ t Av + 


Atf 



and 


+ 9" et Ax nrt + K At AM* + B^Ax^ 


where we can seethe equation of Ax fltt 
page. 


( 22 ) 

at the bottom of the 


□ 


E. System Model 


Based on the above subsystem models of DGl a DG2, ESS 
and network, the small-signal model of a multi -DG micro grid 
system can be obtained. Fig. 6 describes the block diagram 


= A„Ajr + *„Ax- + (17) ii« At'/ Ajjf Ai« At 1 } At} At 1 } Av™ h{ GC f 


Where 
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Fig. 5. Inner current controller. 



Fig. 6. Block diagram of the global small-signal model of a 

multi-DG micro-grid system. 

for developing state-space equations with all of the subsystem 
models integrated together. The node voltages (outputs of 
network) are taken as the input variables to the corresponding 
subsystems, and the node currents (outputs of subsystem) are 
taken as the output variables to the network. 

The complete small-signal model and system state matrix 
(as given in (23)) can be obtained using the subsystem models 
de-noted by (8), (13), (17), (18), (21), and (22): 


&sysAx S yg — H“ -^sysA(j S y S H“ ^sysAtlg^^ 3) 


where 


Ax s , s = [Ax'** Ax” Ax'* A*“ kz net \ 


7 r 


Av syt = [Av w Av c ] 


7 


it All M - 




Au w Au* Au cc 


i7’ 


TABLE I 

System Parameters 


Elements 

Type 

Value 

DG1 

Rotor type 

15 kVA 

DG2 

Rotor type 

31.3 kVA 

ESS 

Inverter type 

16 kVA 

Line 1 

Series RL Branch 

0.05128+/0.008 (fl) 

Line 2 

Series RL Branch 

0.04487+/0.007 (O) 

Line 3 

Series RL Branch 

0.01282+/0.002 (Q) 

Loadl 

Series RL Load 

22.25 kW/18 kvar 

Load2 

Series RL Load 

9.5 kW/1.2 kvar 

Qc P 

Capacitor group 

6 kvar 


TABLE n 

Eigenvalues of the Study Microgrid System 


Eigvj) 

\nlub 

(Ijfaet) 

Im 

(nulfccc) 

vulues 

Real 

(l/sra) 

Im 

(riid/scx;) 

1. 2 

-99} 50 

±38826 

24 

-16006 

6 

L 4 

'W.4H 

±38192 

25 

-IO0OC 

6 

5, 6 

-98.417 

±340.89 

26 

-10060 

0 

7p 8 

-2497.7 

±3U 

27 

-9997 5 

6 

9- 16 

-1992 8 

±114 

28 

-9997J 

6 

11. 12 

-393.32 

±314 

29 

-166.97 

6 

13, 14 

-I6|.l8 

±•105 51 

30 

-60 

6 

35. 16 

-6.57 

±83.250 

31 

-50.831 

6 

17- 18 

-63.241 

±66.15 

32 

-24.409 

6 

19r 20 

■1.3834 

±16.357 

33 

-12.695 

6 

II. 21 

A 2636 

±9,9276 

34 

-2.4 

6 

2? 

-10000 

0 

35 

-6,41829 

6 


III. DYNAMIC STABILITY ANALYSIS OF 

MICROGRID 

The small-signal model can be used for eigenvalues 
anal-ysis to investigate the dynamic stability of microgrid 
system [15]—[18]. System stability margin and corresponding 
operation strategies can be obtained by tracking the loci of the 
eigenvalues with the variation of system parameters or 
steady-state oper-ating state. 

A. Eigen Analysis 

The proposed microgrid system shown in Fig. 1 is taken as 
an example with the main parameters listed in Table I. 

Table II shows the eigenvalues of the microgrid study system 
under autonomous operation mode calculated by the linearized 
state matrix formulated by (23). There are 35 eigenvalues with 
all of the real parts being negative. Eigenvalues 1 to 22 are 
eleven pairs of complex-conjugate eigenvalues which represent 
eleven oscillatory modes. Eigenvalues (5, 6) and (17, 18) 
rep-resent the active and reactive power support by DG2 
respec-tively corresponding to the electromagnetical interactions 
be-tween DG2 and the network. Eigenvalues (13, 14) and (15, 
16) represent the active and reactive power supplied by DG1 
re-spectively corresponding to the electromagnetical interactions 
between DG1 and the network. Eigenvalues (19, 20) and (21, 
22) represent the reactive and active power support by ESS 
re-spectively corresponding to the electromagnetical 
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interactions between ESS and the network. The other five 
oscillatory modes represent the electromagnetical interactions 
among the network. The oscillation frequency and the 
damping of them depend on given operation conditions. 

B. Sensitivity Analysis 

Sensitivity analysis of system eigenvalues can be used to 
in-vestigate the dynamic stability margin of microgrid. It is 
con-ducted with respect to changes in various operation 
conditions. 
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Fig. 7. Loci of eigenvalues (5, 6) and (17, 18) when the capacity of 
DG2 is changed from 30 kVA to 80 kVA. 


Fig. 7 shows the loci of eigenvalues (5, 6) and (17, 18) 
where the capacity of DG2 is changed in condition that ESS is 
not equipped and the other system parameters shown in Table 
I are not changed. It reveals that the capacity increase of DG2 
from 30 kVA to 80 kVA lead to the departure of the conjugate 
pairs of eigenvalues (5, 6) and (17, 18) from the right-hand 
plane to the left- hand plane, which increases the damping of 
the two modes and system stability margin is improved. 

The characteristics analysis of DG2 reveals that the active 
and reactive power output of DG2 increase with its capacity, 
which result in the enhancement of system active/reactive 
power sup-port capacity. This is consistent with the trend of 
eigenvalues (5, 6) and (17, 18) based on sensitivity analysis. 
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Fig. 8. Loci of eigenvalues (5, 6) and (17, 18) when the 
capacity of DG1 is changed from 4 kVA to 20 kVA. 


Fig. 8 shows the loci of eigenvalues (5, 6) and (17, 18) 
when the capacity of DG1 is varied in condition that ESS is 
not equipped and the wind speed is steady at 10 m/s. It can be 
seen that with the capacity of DG1 increasing from 4 kVA to 
20 kVA, eigenvalues (5, 6) move to the left -hand plane, 
which means the damping enhancement of the mode and 


system stability margin improvement. However, eigenvalues 
(17, 18) move to the right-hand plane, resulting the damping 
decrease of the mode and system stability deterioration. 

The sensitivity analysis of DG1 reveals that the active 
power output of DG1 rises with the increase of its capacity 
under the same wind speed status, meanwhile the reactive 
power absorp-tion of DG1 increase too. It results in the fact 
that system active power support capability is enhanced, 
whereas reactive power support capability is weakened. This 
restricts the install capacity of renewable energies. 

Wind speed variation and disturbance will cause the output 
power fluctuations of DG1, which may affect the stability of 
microgrid. Fig. 9 depicts the loci of eigenvalues (5, 6), (13, 14), 
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Fig. 9. Loci of eigenvalues (5, 6), (13, 14), (15, 16) and (17, 18) when wind 
disturbance occurs with the same DG1 capacity. 


(15, 16) and (17, 18) with wind speed variation in the conditions 
that the capacity of DG1 remain unchanged and there is no ESS 
equipped. It indicates that the increasing of wind speed from 7 
m/s to 11 m/s result in the departure of eigenvalues (13, 14) and 
(15, 16) to the right-hand plane, which will reduce the damping 
of the mode and then system dynamic stability. 

It can also be seen from Fig. 9 that the increase of wind 
speed results in the departure of eigenvalues (5, 6) to the 
left-hand plane and eigenvalues (17, 18) to the right-hand 
plane. Ac-cording to the characteristics of squirrel 
asynchronous generator of DG1, the active power output 
increases with the wind speed, while the reactive power 
absorption of DG1 also increases. This results in that system 
active power support capacity is enhanced, whereas the 
reactive power support capacity is weakened. This means that 
microgrid dynamic stability differs with the wind conditions. 

The impact of the output power fluctuating of DG1 on 
micro-grid can be mitigated or eliminated through the equipment 
of ESS. With proper control of power converter between battery 
bank and the microgrid, ESS can operate flexibly, 
absorbing/re-leasing energy quickly to maintain the 
instantaneous power bal-ance of the system. As a result, the 
stability is improved greatly. 

Fig. 10 shows the loci of eigenvalues (5, 6) and (17, 18) when 
wind speed varies from 6 m/s to 11 m/s with ESS equipped in the 
study microgrid system. ESS is controlled to restrain system 
active power fluctuating, and to restrain system active/reactive 
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Fig. 10. Loci of eigenvalues (5, 6) and (17, 18) when wind 
speed varies with ESS equipped, (a) ESS is controlled to 
restrain active power fluctuation, (b) ESS is controlled to 
restrain active/reactive power fluctuation. 


power fluctuating, as shown in Fig. 10(a) and (b), 
respectively. Compared with Fig. 9, the change of eigenvalues 
loci of (5, 6) and (17, 18) in Fig. 10(a) under wind disturbance 
can be re-strained greatly, while as shown in Fig. 10(b), they 
are almost unaffected with wind disturbance. It reveals that 
microgrid dy-namic stability can be improved with proper 
control of ESS. 


IV. SIMULATION AND EXPERIMENT RESULTS 

A simulation platform of the study microgrid system shown 
in Fig. 1 is established under MATLAB/SIMULINK 
environ-ment. Table I shows the main parameters of the 
system. The sim-ulation is implemented in time -domain to 
test and verify system performance under wind gust 
disturbance with/without energy storage. 

The microgrid is controlled to operation in autonomous 
mode. It stays originally in a steady operating state with wind 
speed of 10 m/s. Then a wind gust occurs at t— 340s. Fig. 11 
shows the power output process of DG2 under the gust 
disturbance with/without ESS, respectively. 

As can be seen from Fig. 11(b), the gust disturbance lead to 
DG2 active power fluctuating obviously with the fluctua-tion 
ratio of 10% when ESS is not equipped. However, the ac-tive 
power fluctuating is restrained largely with the fluctuation 
ratio of around 1% and shortened recovery time when ESS is 
used. Meanwhile, the wind gust disturbance lead to DG2 
reac-tive power fluctuating with the fluctuation ratio of 12% 
when ESS is not equipped, and decreased to 1.6% with ESS 
used, as shown in Fig. 11(c). 

The simulation results reveal the fact that ESS can release 
or absorb power quickly and flexibly which track and 
compensate the change of wind power precisely to mitigate 
system instanta-neous power unbalance. ESS absorbs active 
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Fig. 11. System response to gust disturbance (with/without 
ESS), (a) Wind speed, (b) Active power of DG2. (c) Reactive 

power of DG2 


reactive power when wind speed increases, and releases 
active power and absorbs reactive power when wind speed 
decreases. The amount of active/reactive power 
released/absorbed by the ESS can be altered conveniently 
through the control of the VSC between energy storage and 
microgrid. 

A microgrid experimental platform of the same configuration 
and parameters as shown in Fig. 1 and Table I. is established in 
the Microgrid Research and Test Laboratory of IEE, CAS. 

Fig. 12 shows system response to wind gust disturbances 
without ESS. The microgrid system is originally in autonomous 
mode with 6 kW power output of DG1. Following the gust 
oc-curs at £=4QJOb, the output power of DG1 and DG2 change 
accordingly, as shown in Fig. 12(a) and (b). The active power 
output of DG2 decreases when that of DG1 increases with wind 
speed to maintain global active power balance. Meanwhile, the 
reactive power released from DG2 increases too to compensate 
the rising of reactive power absorption of DG1. On the contrary, 
when wind speed decreases, the active power output of DG1 
decreases and that of DG2 increases. Meanwhile, the reactive 
power output of DG2 decreases with the decline of the reactive 
power absorption of DG1. The experiment results show that the 
active/reactive power of DG2 fluctuates obviously during the 
gust disturbance without ESS. This is not favorable for a thermal 
engine since the frequently changing of operation state may 
im-pose extra mechanical and electrical stress on the units. 

Fig. 13 shows system response to gust disturbance with ESS. 
As mentioned above, ESS can track the output power changes 
of DG1 precisely and mitigate the adverse effects through 
active/reactive power output/input of ESS quickly and 
flexibly due to its decoupling controller. As a result, the 
instantaneous power unbalance is mitigated and the frequency 
and voltage quality and stability of the microgrid is improved 
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greatly. Specifi cally, ESS absorbs active power and releases 
reactive power when wind speed increases, whereas output 
active power and absorbs reactive power when wind speed 
decreases, as shown in Fig. 13(a) and (b). According to the 
experiment results, the active/reactive power of DG2 is 
main-tained stable during the gust disturbance with the help of 
ESS. Moreover, the thermal engine in a micro grid can be 
maintained at a high-efficient operating state and system fuel 
economic is improved. 


V. CONCLUSION 

In this paper, a comprehensive small-signal model of a typical 
microgrid is presented including synchronous generator, 
asyn-chronous generator, and voltage source converter with 
corre-sponding control systems. The model covers almost all the 
pos-sible power interfaces of DG and ESS existing in microgrid 
or distributed generation system. All the subsystems are 
modeled individually and then combined on a global reference 
frame. The model of different kind of microgrids can be 
obtained easily based on this paper with proper adjustment of 
the number and parameters of the subsystems. 

The model is analyzed in terms of system eigenvalues and 
their sensitivity to different operating state and control strategies 
of DG or ESS. It is observed that the dynamic stability of an 
au-tonomous microgrid is influenced by system configuration 
and the variability of intermittent generation. It can be improved 
re-markably with ESS equipped due to its ability to mitigate the 


in-stantaneous active/reactive power unbalance. Therefore, 
prop-erly designed and controlled energy storage is an essential 
solu-tion to the reliable operation of microgrid as well the 
effective utilization of renewable energy. 

Results obtained from the model and eigenvalues analysis 
were verified through simulation and experiment on a study 
microgrid system. It can be observed that the model and its 
eigenvalues successfully denote the relationship between the 
dynamic stability and system configuration and the 
corre-sponding operation state. 

This paper provides a theoretical foundation for the 
anal-ysis, confi guration and operation of microgrid 
especially in au-tonomous operation mode. It enables the 
design and control of microgrids to achieve the objective of 
both reliability of power system and effective utilization of 
renewable resources. 
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